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A complete solution of the ground-state problem for an Ising model on the Shastry-Sutherland 
lattice with an additional interaction along the diagonals of "empty" squares in an applied magnetic 
field is presented. A rigorous proof is given that this interaction gives rise to a plateau at one-half 
of the saturation magnetization. Such a fractional plateau has been observed in some rare-earth- 
metal tetraborides, in particular, in strong Ising magnets ErB4 (where it is the only one) and TmB4 
(where it is the broadest one), but its origin has remained unclear. Our study sheds new light on 
the solution of this problem. 

PACS numbers: 75.60.Ej, 05.50.-fq, 75.10.Hk, 75.10.-b 



I. INTRODUCTION 

Here we consider an Ising model on the Shastry- 
Sutherland (SS) lattice with an additional diagonal inter- 
action. The SS lattice is topologically equivalent to the 
Archimedean 3^.4.3.4 lattice (Fig. 1). It has been shown 
that the ground-state problem of the quantum Heisen- 
berg model on the SS lattice in an applied magnetic field 
can be solved exactly^ This model has attracted con- 
siderable interest when the magnetic subsystem of some 
compounds has been shown to consist of weakly cou- 
pled layers of magnetic ions arranged on a lattice that 
is topologically equivalent to the SS one. These com- 
pounds are referred to as Shastry-Sutherland magnets. 
SrCu2 (603)2 is the most known and the most studied 
among them3i2. But during the last decade, many other 
SS magnets have been discovered. In particular, this 
concerns an entire group of rare-earth-metal tetraborides 
KB4 (i? = La - Lu). Magnetic ions R^~^ carry spins large 
enough to be considered as classical Heisenberg ones. If, 
in addition, the crystal field effects are strong, then the 
magnets can be described in terms of the effective spin- 
1/2 model under strong Ising anisotropy. This is the case 
of TmB4 and ErB4. 

SS magnets have attracted particular interest, since 
they exhibit sequences of fractional magnetization 
plateaus. A single plateau at one-half of the saturation 
magnetization {m/ms — 1/2) has been observed in ErB4. 
TmB4, in addition to an extended 1/2-plateau, exhibits 
a sequence of narrow plateaus with fractional values of 
m/ms = 1/6, 1/7 up to 1/12 for temperatures below 4 K 
(the magnetic field is normal to SS planes) 

Some attempts have been made to explain the ori- 
gin of the fractional magnetization plateaus in terms of 
the Ising model on the SS lattice, however, this model 
was shown, both numerically^^?, and analytically,^ to pre- 
dict a single fractional plateau — at m/m^ = 1/3. This 
suggested that only quantum models and maybe with 
longer-range interactions could describe the 1/2-plateau 
and other fractional plateaus observed in the rare-earth- 
metal tetraborides/'^" In Refs.fiol andlTTI. a 1/2-plateau 
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FIG. 1: (a) Archimedean 3^.4.3.4 lattice, (b) Shastry- 
Sutherland lattice and (c) the lattice formed by magnetic 
Cu^^ ions in SrCu2(B03)2. All the three lattices are topolog- 
ically equivalent, (d) An extended Shastry-Sutherland lattice 
considered here. 



was obtained using the quantum SS model (and also its 
Ising limit) with additional interactions J3 (along the di- 
agonals of "empty" squares) and J4 (the next-nearest 
neighbors along edges). 

In our previous work, we showed that the fractional 
plateaus in an Ising-type model on the SS lattice can 
be generated by the long-range RKKY-type interactions 
(these are present in rare-earth-metal tetraborides since 
the latter are good metals), and that the 1/2-plateau 
is given rise by the additional diagonal interaction 
Recently, a 1/2-plateau (as well as a narrow 2/3-plateau) 
was shown numerically to be generated by dipole-dipole 
interactions 

In our opinion, to explain the sequence of the mag- 
netization plateaus in TmB4, it is sufficient to find the 
ground states of an Ising-type model on an extended SS 
lattice with long-range interactions. It is difficult to find 
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a numerical solution but our analytical method for deter- 
mining the ground states of lattice gas models or equiv- 
alent spin ones makes it possible. 

Here, we present a complete solution of the ground- 
state problem for the Ising model on the SS lattice with 
additional interaction J3 along the diagonals of "empty" 
squares (in what follows, we refer to this model as the 
Ising model on the extended SS lattice). To find this so- 
lution, we use the method of basic rays and basic sets of 
cluster configurations which was developed in our previ- 
ous works. Moreover, we generalize the method, and 
consider configurations of two (instead of one) different 
clusters that is quite natural for the SS lattice with ad- 
ditional diagonal bonds. 

We rigorously prove the existence of a rather wide 1/2- 
plateau in this model. It corresponds to three different 
phases and appears under the condition that the interac- 
tion Ji along the edges of squares and at least one of the 
two remaining interactions — along the SS diagonals (J2) 
and along the diagonals of the "empty" squares ( J3) — is 
antiferromagnetic. It depends on the signs of the interac- 
tions J2 and J3 which of these three phases is realized. A 
single 1/2-plateau exists, in particular, in the case when 
the interactions Ji and J2 are antiferromagnetic and J3 
is ferromagnetic (under the condition J2 < — 2J3). A sin- 
gle 1/2-plateau has been observed in ErB4, but its origin 
remains unclear. Hence, our investigation sheds light on 
this problem. 

The paper is organized as follows. In Sec. II, a two- 
cluster approach to the determination of the ground 
states of Ising-type models is developed and a solution 
of the Ising model on the extended SS lattice is found. 
In Sec. HI, the full-dimensional structures are deter- 
mined and analyzed. In Sec. IV, the structures at the 
three-dimensional boundaries of the full-dimensional re- 
gions are considered and applied to predict the order of 
phase transitions between full-dimensional phases. These 
are also important for the investigation of the effects of 
longer-range interactions. In Sec. V, the ground states at 
the two-dimensional boundaries of the full-dimensional 
regions are analyzed. In Sec. VI, the ground-state phase 
diagrams in the (/i, J2)-plane are presented as well as pos- 
sible sequences of phases for the magnetization processes 
in ErB4 and TmB4. Section VII gives some conclusions. 



II. SOLUTION OF THE GROUND-STATE 
PROBLEM FOR THE ISING MODEL ON THE 
EXTENDED SS LATTICE: A TWO-CLUSTER 
APPROACH 

In our previous studies )^ '^'^'^^ we have constructed the 
ground-state structures for Ising-type models (or equiva- 
lent lattice-gas models) with configurations of some clus- 
ter. We used only one kind of cluster. For instance, a 
cluster in the form of a triangle with an SS diagonal as the 
hypotenuse is sufficient for the Ising model on the con- 
ventional SS lattice;S Here we consider the same model 




FIG. 2: (a) Clusters for the determination of the ground 
states: (left) a cluster with an SS diagonal (interaction J2) 
and (right) a cluster without SS diagonal, but with ordinary 
ones (interaction J3). (b) The energy distribution of a site 
between four squares sharing this site. 



but with additional interaction J3 along all the diagonals 
of the squares without SS bonds [Fig. 1(d)]. To construct 
the ground-state structures for this model one small clus- 
ter is insufficient; one should use two different clusters: 
a square with an SS diagonal and a square without it 
[Fig. 2(a)]. However, there is no crucial difference in com- 
parison with the one-cluster approach. The main idea is 
the same and consists in introducing "free" coefficients 
which account for the fact that the energy contribution of 
sites and some bonds can be distributed between clusters 
in various ways. (As will be clear later, "free" coefficients 
make it possible to equalize and minimize the energies of 
two or more cluster configurations at once). Each edge 
belongs to two squares (one of each type) and each site 
belongs to four squares (two of each type). The energy 
contribution of an edge (site) can be distributed in var- 
ious ways between the squares sharing this edge (site). 
The way of distribution is determined by the "free" co- 
efficients a, /3, and 7 in the expressions for the energy 
contributions of different squares [see Fig. 2(b)] 

e — -I- aJi{ai(j2 + cr20'3 + o'3CT4 + a'4cri) 

- -lh[P{cJl + 0-3) + (1 - ^)((T2 + 0-4)], 

e = J3(crio-3 + 0-20-4) + (1 - a) Ji(ctio-2 + 0-20-3 

1 - 7 

+ O-3O-4 + O-4O-1) —h{(Ti + 0-2 + 0-3 + 0-4).(l) 

Here e and e are the energy contributions of the squares 
with SS bond and without SS bond, respectively; Ji, 
J2, and J3 are the parameters of interaction along the 
edges, SS diagonals, and ordinary diagonals (only for 
the squares without SS diagonals), respectively; h is the 
applied magnetic field. The numbers ai, 0-2, 0-3, and 
0-4 {cTi = ±1) define the spin configuration of a square. 
There are nine configurations of a square with an SS di- 
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TABLE I: Basic rays and basic sets of configurations for the Ising model on the extended Shastry-Sutherland lattice. 





Basic ray 


Basic set 


Full-dimensional 




"Free" 






{h, Ji, J2, J3) 


of configurations 


structures 




coefficients 




ri 


(0,0,0,-1) 


00 o# ## 
00 mo •• 


1, I, 3, 9, 10, To 




Arbitrary 




r2 


(0,0,-1,0) 


00 AO OA AO OA AA 

ao do mo ai mi mi 


1, 1, 3, 7, 8 




Arbitrary 




ra 


(0,-1,0,1) 


00 #0 

00 AO 


1, 1, 5, 8 




a = 




r4 


(0,-1,2,0) 


00 00 AO AA AA 

VJAJ I/-* "J* 

do ¥0 «fo d« 


1, 1, 5, 9 




a = 1 




rs 


(0,1,2,0) 


00 AO OA AO AA 

vj—' '-'J* 

*o *o ifo cJ« d« 


3, 4, 4, 5, 9 




Q = 1 




re 


(2,0,0,1) 


#0 ## ## 

•0 o« •• 


1, 5, 6, 7, 8 




7 = 




rr 


(4,1,0,0) 


OA AO AA OA AA 


1, 3, 4, 6, 7, 10 


a 


= 1, /3=i7 = 


1 


rs 


(1,0,1,0) 


00 AO OA AA OA AA 


1, 4, 5, 6, 9, 10 




/3 = 1,7=1 




rg 


(2,1,0,1) 


#0 o# ## 

• AO O* 


3, 4, 5, 6, 7, 8 




Q = 0, 7 = 






f2 1 2 ~1) 


mpqmmp mm^f q»^= i, 00 om mm ••+0 
• mo om om mm w 00 mo mm om 


3, 4, 9, 10 


(X 
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re 


(-2,0,0,1) 


00 00 mo 
00 mo mo 


1, 5, 6, 7, 8 




7 = 






(-4,1,0,0) 


op op ma om mo 

00 mo ao mo ai 


1, 3, 4, 6, 7, 10 


a 


= 1, /3 = |,7 = 


1 


rs 


(-1,0,1,0) 


qpopmpmp mp mm 

00 mo ao mo ai ai 


I, 4, 5, 6, 9, To 




/? = 1,7 = 1 






(-2,1,0,1) 


mo om 00 
mo mo mo 


3, 4, 5, 6, 7, 8 




a = 0, 7 = 






(-2,1,2,-1) 


mpqmmp qp^ mp^, ,, po om mm ao^h 
«p «p cf» «p Qsa ||pp«p»»»P 


3, 4, 9, To 


a 


= 1, /3=i,7 = 


1 



"Configurations marked by asterisk enter the structures in blocks shown in Fig. 4. 
''Condition symmetric to the condition for rig. 



agonal and six configurations of a square without SS di- 
agonal: §0, *§, «§, *?, S§, dZ, S5, iZ II PP, So, Sp, So, 

oSj [open and solid circles denote spin down (—1) and 
spin up (+1), respectively; the squares of different types 
are separated by the symbol ||]. 

As in our previous studies, here we use the method of 
basic rays (vectors) and basic sets of cluster configura- 
tions. In brief it is as follows. 

Consider the parameter space (ft., Ji, J2, J3) of the 
model. Each ground-state structure corresponds to a re- 
gion in this space. The region can be presented as a 
solution of a set of uniform linear inequalities; therefore, 
it is a polyhedral cone. Let us recall that a polyhedral 
cone is a conical hull — that is, all the linear combina- 
tions with nonnegative coefficients — of a set of vectors. 
It is completely determined by its edges or vectors along 
them. Among the polyhedral cones corresponding to all 
possible ground-state structures of the model, the most 
important are those with the dimensionality equal to the 
dimensionality of the parameter space. We refer to these 
regions and corresponding structures as full-dimensional. 

If a structure is a ground-state one at two points of the 
parameter space, then it is a ground-state structure in the 
entire line segment connecting these points. This is a con- 
vexity property; the region corresponding to a structure 
is always convex. The convexity property makes it possi- 
ble to determine the ground-state structures in any point 
of the parameter space if all edges of the full-dimensional 
regions (basic rays) and all ground-state structures in 
these edges are known. 

First we determine the cluster configurations which 
generate ground-state structures and then the structures 



themselves. By definition, a set of cluster configura- 
tions generates a structure if all the cluster configura- 
tions in the structure belong to the set. A set of cluster 
configurations generates all the structures at the point 
{h, Ji, J2, J3) of the parameter space if the following three 
conditions are satisfied: (1) values of the "free" coeffi- 
cients a, (3, and 7 in Eqs. (HJ can be chosen in such a 
manner that all the configurations of the set have the 
same energy, (2) this energy is lower than the energies 
of all the remaining configurations, and (3) at least one 
structure can be generated by the configurations of the 
set. It should be noticed that this group of conditions is 
sufficient but not necessary; only the third condition is 
necessary. 

The basic rays and the corresponding sets of cluster 
configurations (basic sets) of the model under consider- 
ation are listed in Table I. Configurations of clusters of 
different types are separated by the symbol ||. If config- 
urations of a type of cluster are not indicated in a set, 
then all the configurations of this type belong to the set 
(i.e., the configuration of this type of cluster can be ar- 
bitrary) . The last column of Table I contains the values 
of the "free" coefficients for which the configurations be- 
longing to the basic set have the same energy which is 
lower than the energies of the remaining configurations 
(in the corresponding ray, of course) . 

Having determined the full-dimensional regions and 
structures on the basis of Table I, we can show that the 
set of basic rays is complete. This will be done later. It 
should be noticed that we present here only the final re- 
sults, the correctness of which can be easily verified. In 
reality, first we found the full-dimensional ground-state 



structures of the model (all these but one are determined 
in Ref. ^ and only then we calculated the basic rays. 

In all the basic rays, except for rio and r^Q, all the 
structures generated by the configurations of the corre- 
sponding basic set are ground-state structures in these 
rays. As to the rays rio and r^Q, an additional condition 
should be satisfied. For the ray rio this condition can be 
formulated as follows: the configurations marked by the 
asterisk enter the ground-state structures in the blocks 
shown in Fig. 4 and these blocks do not overlap either 
with themselves or with other squares. (For the symmet- 
ric ray r]~o the symmetric condition should be satisfied). 
Let us prove this. 

Using Eq. ([Ij, one can calculate the energies of all 
the configurations of the clusters at the point h — 2Ji, 
J2 = 2Ji, and J3 — — Ji (the ray rio) for the following 
values of the "free" coefficients: a = 1, (3 = 1/2, and 
7 = 1. 

S: (-2Ji); SS, S§ (O); S§, S: (2Ji); K (4Ji); S§ 

(6J1); g§ (12Ji); 88, 2S, (0); 28, S: (2Ji); :8 (4Ji). 
The zero energy level is shifted here by 2Ji > 0. The 
energy value indicated after each group of configurations 
is the same for all the members of the group. We refer to 
the squares with negative, zero, and positive energies as 
"negative," "zero," and "positive" squares, respectively. 

Let us prove that, in the ray rio, the energy of an ar- 
bitrary structure cannot be negative. Near each negative 
square 55, in the region shown in Fig. 3(a), at least one 
square among the three squares of the region is positive. 
Let us group all the squares in a structure in the follow- 
ing way: each zero square forms a group by itself; each 
negative square enters the group of the positive square 
situated near it in the region shown in Fig. 3(a). If a 
negative square enters two or three groups, then its en- 
ergy is distributed into equal parts between these groups. 
It it easy to see that, for such a grouping rule, positive 
squares §§, and 28 cannot form a group with any neg- 
ative square. Squares §§ and «8 can form groups with 
one negative square only, but the energy of the groups is 
positive. Square H can form a zero-energy group with 
two negative squares. In this case, however, other groups 
with positive energies are formed inevitably [Fig. 3(b)]. 
The square 2* can form a sigle zero-energy group with a 
negative square in the way shown in Fig. 3(b). The only 
square that can form a group with negative energy (— Ji) 
is the square *• [Fig. 3(d)]. However, then there appears 
one of the two positive squares, 58 or *•■ The first one 
forms a group with the energy 3 Ji which cannot be com- 
pensated by the negative energy of the group formed by 
the square *•• The second square *• can form a zero- 
energy group, but in this case the situation is repeated 
and so the negative-energy group of the first square ** 
generates an infinite half-stripe in which all the remain- 
ing groups have zero energies. Hence, the number of 
negative-energy groups in the structure [Fig. 3(d)] can 
be infinitesimal only. To summarize, the ground-state 
structures in the ray rio, can contain, along with zero 
squares, the positive squares 25 and *• but only in com- 
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FIG. 3: (a) The region (solid line) near a negative square 
containing at least one positive square, (b) Square can 
form a zero-energy group but some positive-energy groups 
emerge then inevitably, (c) Zero-energy group with a square 
2*. (d) Square *• can form a negative-energy group (solid 
line), but in this case a new square *• (dotted line) emerges 
with a group that should have zero energy and, therefore, this 
square yields the same square and so on. 




FIG. 4: Blocks with the squares 2« and *• for the ground- 
state structures in the ray no. 



binations with the square §5 shown in Fig. 4. 



III. FULL-DIMENSIONAL GROUND-STATE 
STRUCTURES AND REGIONS 

Table I represents a complete solution of the ground- 
state problem under consideration. Using this Table and 
the convexity property, it is easy to determine the full- 
dimensional ground-state structures and regions. All one 
has to do is to find the subsets of the basic sets of config- 
urations with the following properties: (1) each of them 
is a subset of at least four basic sets and the linear hull of 
the corresponding basic vectors is full-(four-)dimensional; 
(2) the cluster configurations of such a subset generate 
at least one structure. 

The full-dimensional regions and the corresponding 
subsets of configurations which generate all the ground- 
state structures in these regions are listed in Table II. 
The first column gives the numbers of full-dimensional 
regions. Taking into account the symmetry of the model 
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TABLE II: Full-dimensional regions and ground-state structures of the Ising model on the extended Shastry-Sutherland lattice 
(for /i > 0). 



Region 


Generating 




Characteristics 


Magneti- 






Number 


Condition for existence 


(Structure) 


configurations 




of the structure 


zation 




Basic rays 


of 3-faces 


in the plane {h, J2) 


1 


K II 22 
•• II •• 




2 /, + ^ + 7-, - h 

[1 II 1] 


1 


r7, 


1*0 T/I Vi^ 
^2;A37^45^D; 


8 


Always 


r> 
o 


•o a» II 90 




9 7"-, -1- ■^^ 1 r 

[illl 1] 


n 
u 


fi, 
riQ 


"2, ^5, r7, rg, 

1, 1*7 , , r-]^Q 


1 1 

i. J. 


Jl u, J3 ^ Ji 


4 


0» M II 0« 
• O 0« II vo o* 




3 6 3 3 

ri 2 II 1 2-1 

L3 ' 3 II 3 ' 3 J 


1 

i/O 




TTi-r Vc 1*1-1 1*1/-! 


5 


IJ3I < Jl 





56 II 2g 
•o II mo 




"'2 7, 
2 3 

[1 II 1] 


u 


1'3, 


"4, r5, re, rg, rg, 
' ^8 ' ^9 


q 


J3 ^ U 


g 


•o o* 1 




[o,iio|| 1] 


1 /9 


re, 


17, 185 J- 9 


A 




7 


SSS222S2 II 


om 


J2 ft 
2 2 

[Desorder |[ 1] 


1 /2 


''2, 


rs ; 1*7, rg 
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7i > n 7o < n 7q > n 


8 


§§SSS§Kll 


•o 
•o 







r2, 


r3,r6,r9,r^,rg 


6 


J2 sS 0, J3 ^ |Ji| 








L4 ' 4 ' 4 ' 4 11 ^1 
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•P 11 oo o« 
mo II oo mo 




-# + Ja 

fl II i i il 

L II 4 ' 2 ' 4l 





ri, 


r4,r5,r8,rio, 
I ''10 


8 


J2 ^ 2|Ji|, J3 < 


10 


•» o^ II o« 

Cf* II AO 






1/2 


ri, 


r7,rg,rio 


4 


J2 > 0, J3 s; 



i II i i] 

L2 ' 2 II 2 ' 2J 



with respect to the field inversion with simultaneous flip 
of all spins, we indicate only the regions (structures) 
with zero and positive magnetization. We denote a re- 
gion (structure) with negative magnetization by the same 
number as the symmetric region (structure) but with a 
bar over the number. The third column of Table II gives 
energies per site of structures and, in the square brack- 
ets, the fractional contents of cluster configurations in 
the structures (see Ref. [isl ). If the fractional contents 
of configurations in a structure can vary (region 7), then 
there is a degeneracy marked by the word "Disorder." (In 
region 6, there is degeneracy as well, but the fractional 
contents of configurations do not vary). The fourth, fifth, 
sixth, and seventh columns indicate, respectively: mag- 
netization of the structure, basic rays that define the cor- 
responding region, the number of three-dimensional faces 
of the region, and the conditions for the existence of the 
region in the (ft., J2)-plane. 

The full-dimensional ground-state structures of the 
model under consideration are depicted in Figs. 5 and 
6. They are constructed with the cluster configurations 
listed in Table II. All the full-dimensional structures, ex- 
cept for structures 6 and 7, are fully determinate. In 
the regions 6 an 7, an infinite number of the ground- 
state structures occur: phases 6 and 7 are disordered 



(Fig. 6). It is easy to see that the disorder of phase 7 is 
one-dimensional: it is ordered in one direction and disor- 
dered in another. 

It is more difficult to determine the character (i.e., the 
dimensionality) of disorder in phase 6. In structures of 
phase 6 (and even at the boundary between phases 4 
and 6), the structural element shown in Fig. 7 generates 
an infinite half-stripe. On the basis of this fact one can 
prove that the disorder of phase 6 is one-dimensional. 
The number of squares and S{ is infinitezimal in com- 
parison with the number of remaining squares. Although 
their energies are equal to the energies of other squares 
which generate structure 6, these squares represent a kind 
of zero-energy defects. Hence, structure 6 can be consid- 
ered as a simple mixture of structures 6a and 6& depicted 
in Fig. 5. 

It is interesting that the full-dimensional structure 9 
is chiral. This is better seen on the lattice shown in 
Fig. 1(c). Two different structures 9 are possible: a struc- 
ture with the SS bonds twisted counterclockwise around 
"open" and "solid" squares and its chiral twin with the SS 
bonds twisted clockwise (Fig. 8). Hence, the interaction 
J3 lifts the degeneracy of the Ising dimer phase (opposite 
spins on the SS diagonals)^ and thus two ordered phases 
are produced; one of these is chiral. It should be noticed 
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FIG. 5: Full-dimensional ground-state structures of the Ising model on the extended SS lattice (for h 0). Phases 1, 3, 4, 
and 5 are, respectively, the fully polarized phase, the Neel phase, the 1/3-plateau phase (or the UUD phase), and the coUinear 
phase. Phases 8 and 9 are chessboard phases. Phases 10, 6a, and 66 are 1/2-plateau phases. A simple mixture of structures 6a 
and 66 is the ground-state structure for region 6. Unit cells are indicated. 





FIG. 6: Full-dimensional ground-state structures of the Ising 
In structure 6 the squares- "defects" are shown. 



model on the extended SS lattice. Disordered phases 6 and 7. 



that this structure is not chiral on the Archimedean lat- 
tice [Fig. 1(a)] with equal interactions Ji and J2. Chi- 
ral structures are ubiquitous in nature, however, their 
emergence as a result of spontaneous symmetry breaking 
remains unclear 



IV. GROUND-STATE STRUCTURES IN THE 
3-FACES OF THE FULL-DIMENSIONAL 
REGIONS 

Let us consider now the ground-state structures in the 
three-dimensional faces (3-faces) of the full-dimensional 
regions. A 3-face is defined by a subset of the set of basic 
vectors for a full-dimensional region. An algorithm for 
the determination of 3-faces of a full-dimensional region is 
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FIG. 7: A structural element in phases 4 and 6 (as well as in 
their boundary) and an infinite half-stripe generated by this 
element. 





FIG. 8: Chirality of structure 9. The SS bonds twist counter- 
clockwise around "open" and "solid" squares (left structure) 
or clockwise (right structure) but in the same manner for all 
the squares. 



described in Ref . Il3l . It is a simple mathematical problem 
and we do not reproduce the description of this algorithm 
here. Having found the 3-faces of all the full-dimensional 
regions, one can easily prove the completeness of the set 
of basic vectors. It is sufficient to prove that each 3-face is 
a 3-face for two full-dimensional regions and, hence, the 
full-dimensional regions fill the whole parameter space 
without gaps and overlaps. 

Knowing the basic vectors for a 3-face, it is easy to de- 
termine the set of cluster configurations which generates 
all the ground-state structures in this 3-face. This set is 
the intersection of the basic sets of configurations for the 
basic vectors of the 3-face. Then one can describe the 
corresponding ground-state structures. 

Basic vectors for the 3-faces of full-dimensional regions 
as well as corresponding sets of cluster configurations 
are presented in Table III. In this Table, a 3-face is de- 
noted by the numbers of full-dimensional regions (phases) 
which share this 3-face. If a configuration of a set is not 
compatible with other configurations of the set (i.e., the 
configuration cannot enters any ground-state structure 
in this 3-face), then it is separated by the symbol — . 
The character of transition between two regions (phases) 
depends on the ground-state structures at their bound- 
ary (i.e., at their common 3-face). If the set of configu- 
rations for a 3-face generates only the structures which 
are the ground-state ones for the full-dimensional regions 
sharing this 3-face and nothing else, then there is a first- 
order phase transition between these regions (notation 
"Jump" in Table III). This kind of phase transitions ex- 
ists between phases 1 and 1, 1 and 3, 1 and 9, 3 and 
9, as well as 4 and 10 (and also between pairs of sym- 
metric phases). If for any value of magnetization from 
the interval between the values of magnetization for two 



neighboring full-dimensional phases, at least one ground- 
state structure can be constructed, then the transition 
between these phases is continuous (notation "Cont." in 
Table III). As it is clear from Table III, the most part of 
transitions just possess this property. The situation be- 
tween phases 4 and 9 is more complicated: there are both 
discontinuous and continuous phase transitions between 
them. 

An interesting and important question concerns 
the disorder and entropy at the boundaries of full- 
dimensional regions, in particular, in the 3-faces. Let 
us consider, for instance, a 3-face between the regions 1 
and 10. A typical example of a structure in this 3-face 
is shown in Fig. 9(a): the structure consists of diagonal 
spin-down chains (along J3 bonds) separated by diagonal 
spin-up chains whose number can be arbitrary odd but 
not less than three. This is a simple mixture of structures 
1 and 10. It is clear that order exists along these chains 
but there is disorder in the perpendicular direction. Such 
a disorder can be called one-dimensional. It does not lead 
to macroscopic degeneracy: the entropy per site tends to 
zero if dimensions of the lattice tend to infinity. A sim- 
ilar simple mixtures of full-dimensional structures exist 
also at boundaries of phases 1 and 5, 1 and 8, 3 and 
4, 3 and 8, as well as 3 and 10. At the boundary be- 
tween phases 9 and 10, the ground-state structures are 
not a simple mixture of full-dimensional structures 9 and 
10 but rather a kind of their hybrid, although the disor- 
der is one-dimensional there. They look like structure 10 
with additional spin-down diagonal chains normal to the 
similar chains of structure 10 [Fig. 9(b)]. 

Another kind of disorder occurs, for instance, at the 
boundary between phases 1 and 6. The ground-state 
structures at this boundary can be obtained from a mix- 
ture of structures 6a and 66 (which is a ground state in 
region 6) by flipping a part of spins down [Fig. 9(c)]. It is 
clear that this disorder is two-dimensional since infinite 
number of local changes can be made in the structures 
within the bounds of the ground state. Every fifth spin 
in the structure depicted in Fig. 9(c) is "free," that is, 
can be directed downward as well as upward. Hence, the 
entropy per site is equal to |ln2 or maybe more, since 
Fig. 9(b) does not exhaust all structures. 

In Table HI, the dimensionality of disorder is indicated 
in the fourth column. To find it is sometimes not so 
easy as in the cases considered above. For instance, we 
could not determine the dimensionality of disorder at the 
boundary between phases 4 and 5 (see Fig. 10). The 
disorder in this 3-face is most likely one-dimensional but 
it still should be proved. Let us analyze the disorder at 
some other 3-faces where such analysis is nontrivial. 

One can see from Fig. 11 that at the boundary between 
phases 3 and 5 as well as 5 and 8 the disorder exists in 
two perpendicular directions independently. However, at 
the first boundary the disorder is two-dimensional (and 
the entropy per site is nonzero) and at the second one 
the disorder is one-dimensional. The structures at the 
boundary between phases 5 and 8 are interesting in the 
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TABLE III: Basic rays and ground-state configurations for the 3-faces of the fuU-dimensional regions for the Ising model on 
the extended SS lattice {h ^ 0). The dimensionality of disorder is indicated for continuous transitions. 



Basic rays 


Ground-state configurations 


Transition 


Conditions for existence 


Regions of the 3-face 


for the 3-face 


between phases 


in the plane {h, J2) 



1, 


1 


ri 


r2 


rs, r4 


1, 


3 


ri 


r2 


r? 


1, 


5 


rs 


r4 


rg.re 


1, 


6 


re 


r? 


rg 


1, 


7 


r2 


re 


r? 


1, 


8 


r2 


rs 


re 


1, 


9 


ri 


r4 


rg 


1, 


10 


ri 


r? 


rg 


3, 


4 


rs 


rg 


r7,rio 


3, 


5 


rs 


rg 


^9 


3, 


7 


r2 


r? 


rg 


3 


8 


1 2 


i 9 


^9 


3, 


9 


ri 


ri( 


.r5,r-, 


3, 


10 


ri 


rr 


rio 


4, 


5 


rs 


rg 


r9 


4, 


6 


ry 


rs 


rg 


4, 


9 


rs 


rg 


rio 


4, 


10 


r? 


rg 


rio 


5, 


6 


rs 


rg 


rg 


5, 


8 


rs 


re 


r9,rg , 


5, 


9 


r4 


rg 


rs.rg 


6, 


7 


re 


r? 


rg 



9, 10 



00 11 00 
00 *• II 00 

mo ai mi \ mi W •o \ 00 
mp mm , op mm n mo mm 
«o «• I ¥0 d* II vo 
om mm om mm n 1 mo 

•o «• «• II o* I vo 
•o 
•o 
•o 
•o 



om mp om mm mm mm 

mo ai mi mi II o* 

•p I op ^ II 00 o* 



r2,re,r9 
ri,rg,rio 



00 mo mm 

om mm , 00 
mo mm I 00 



om mp mm n om mm 

mo aiai II •© o» 



I Jo 0§ S» 



•o om 
mo mo 

om mp om mm n o* 1 mo 

«0 «• II MO 0« I VO 

OD mp om mp om mm n mo om 

do do mo ai mi mi II •© #0 
mp om mp n 00 o* 
«ro i(o c<« II 00 00 

o* I 00 
•o I 00 



op mp om mm 
•o *o •o d» 

II 



II 90 0« 

II 90 90 om 
II o* I mo 
II mo om I mo 
00 om mm mm^a 

00 mo mm om 
om mm mm* , 00 

1 mo mm om I 00 



•o 



•o 

•o o* 



op mp am mp am mm 

do do mo di mi mi 



Jump 

Jump 
Cont. (1) 
Cont. (2) 
Cont. (2) 
Cont. (1) 

Jump 
Cont. (1) 

Cont. (1) 

Cont. (2) 
Cont. (2) 
Cont. (1) 

Jump 
Cont. (1) 

Cont. (?) 
Cont. (1) 
Jump -j- Cont. (1) 
Jump 



Ji < 0, J2 > min[-2Ji, -2(Ji + Js)], Js < | Ji| 
Ji > 0, J2 < 0, Js < 
Ji < 0, J2 > max[0, -2(Ji + Js)], J3 > 
Ji > 0, J2 > 0, Js > 
Ji > 0, J2 > 0, Js > 
< - Ji < Js, J2 < 
Ji < 0, J2 > -2Ji, Js < 
Ji > 0, J2 > 0, Js < 

max(0, -2Js) < J2 < min[2Ji,2(Ji - Js)], 
Ji > 0, I Js| < Ji 
J2 = 2(Ji - Js), < Js < Ji 
J2 < 0, < Js < Ji 
J2 < 0, Js = Ji 
Ji > 0, J2 > 2Ji, Js < 
Ji > 0, < J2 < min(2Ji, -2J3), Js < 

J2 > 2(Ji - Js), < Js < Ji 
J2 > 0, < Js < Ji 
J2 > 2Ji, -Ji < Js < 
J2 > -2Js, -Ji < Js < 



I • 



•o 

•o o* 



00 mo I 



Cont. 


(2) 


< Ji < Js, J2 > 


Cont. 


(1) 


J2 = 0, Js > |Ji| 


Cont. 


(2) 


J2 > 2|Ji|, Js = 


Cont. 


(2) 


Ji > 0, J2 = 0, Js > 


Cont. 


(2) 


< Ji < Js, J2 < 


Cont. 


(1) 


Js < -Ji < 0, J2 > 2Ji 



"Configurations marked by asterisk enter the structures in blocks shown in Fig. 4. 



FIG. 9: Examples of the ground-state structures at the boundaries between phases (a) 1 and 10 as well as (b) 9 and 10. The 
disorder at theses boundaries is one-dimensional, (c) An example of ground-state structures at the boundary between phases 
1 and 6. Each gray circle can be either open or solid, hence, the disorder is two-dimensional at this boundary. If all the gray 
circles are open, then we obtain a mixture of structures 6a and 66, which is a ground-state structure in region 6. 
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FIG. 10: An example of a ground-state structure at the 
boundary between phases 4 and 5. It seems that the disorder 
is one-dimensional, however we cannot prove it as yet. 



sense that they are determined by a single condition: ah 
the "empty" squares should have configuration Jo- 

An interesting set of structures also occurs at the 
boundary between phases 4 and 9. The disorder at this 
boundary is one-dimensional, since the structure is com- 
pletely determined by a zigzag-stripe an example of vifhich 
is depicted in Fig. 12. The structure closest to structure 
9 in this set is an ordered structure with m/ms = 1/5 
[Fig. 13(a)]. Hence, there is a jump between structure 
9 and this structure and then a continuous transition 
to phase 4, that is, at the boundary between phases 4 
and 9 there occur both a jump and a continuous tran- 
sition. In addition to the structures of the type shown 
in Fig. 12, there is one structure of another type at the 
boundary between phases 4 and 9. This is a structure 
with m/ms = 1/4 shown in Fig. 13(b). All these struc- 
tures (except for structure 4) are chiral (see Fig. 8 and 
its explanation in the text). 

Another set of structures, where it is difhcult to deter- 
mine the dimensionality of disorder, occurs at the bound- 
ary between phases 4 and 6. An example of a struc- 
ture at this boundary is shown in Fig. 14. It seems that 
the disorder is two-dimensional, however, a more pro- 
found analysis shows that the disorder is nevertheless 
one-dimensional since, as in phase 6, the structural el- 
ement shown in Fig. 7 generates an infinite half-stripe. 

At the boundary between phases 3 and 4 a collection 
of stripe structures exists. Most probably, some of these 
give rise to fractional magnetization plateaus in TmB4. 
We do not describe these structures here, since this is 
done in Ref. The interaction J3 does not lift the de- 
generacy at this boundary. 



'X >■( -X <X s^C si^X X "A v'C 

ag)i(p)irp^ 




FIG. 11: Examples of ground-state structures at the bound- 
ary between phases (a) 3 and 5 as well as (b) 5 and 8. 
Macroscopic degeneracy (residual entropy) occurs only at the 
boundary between phases 3 and 5 but not at the boundary be- 
tween phases 5 and 8 although in the latter case the disorder 
exists in two perpendicular directions independently. Some 
sites are connected with lines for better visualization. 



V. GROUND-STATE STRUCTURES AT THE 
2-FACES OF FULL-DIMENSIONAL REGIONS 

Although the ground states in the 2-faces of full- 
dimensional regions are not so important as the ground 
states in the 3-faces, let us consider them for complete- 
ness. We denote a 2-face by two basic vectors (which gen- 
erate it) in curly brackets. In Table III, vectors of each 3- 
face are ordered in the way that each pair of neighboring 
vectors (the first and the last are neighbors as well) gener- 
ates a 2-face. The set of ground-state configurations for a 
2-face is an intersection of the basic sets of configurations 
for these vectors. The 2-faces of full-dimensional regions 
for /i ^ and the ground-state configurations for them 
are presented in Table IV. The third, fourth, fifth, and 
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FIG. 12: An example of ground-state structures at the bound- 
ary between phases 4 and 9. The rectangles formed by 
the ferromagnetic chains together with "open" and "solid" 
squares are shown as well as one zigzag-stripe that generates 
the whole structure. The structure is chiral: the SS bonds 
twist clockwise around "open" squares and counterclockwise 
around "solid" ones (and vice versa for the chiral twin). 



sixth columns in Table IV indicate, respectively, the di- 
mensionality of disorder in a 2-face, the full-dimensional 
regions sharing this 2-face, the coordinates of the point 
that corresponds to this 2-face in the plane (/i, J2), and 
the conditions for the existence of such a point in this 
plane (if the coordinates are not indicated in parenthe- 
sis, then the 2-face does not intersect the plane (/i, J2) or 
lies completely in it for some values of Ji and J3). 

It is easy to determine the dimensionality of disorder 
in the most part of 2-faces, taking into account that the 
disorder in a 2-face cannot be lower than the disorder in a 
3-face bounded by this 2-face. This task is difficult, how- 
ever, for some 2-faces. Let us consider such 2-faces. In 
the 2-face {r3,r4}, the disorder is two-dimensional. This 
is clear from Fig. 15. It is more difficult to show the disor- 
der in the 2-face {r5, rio} is two-dimensional too. This is 
illustrated in Fig. 16, where the zigzag-strips going from 
left to right can both descend and ascend a little. In the 
2-face {r8,rio} the disorder is two-dimensional. This is 
clear from Fig. 17, where each pair of gray circles should 
contain one open and one solid circle. 

Finally, let us note that the disorder is two-dimensional 
in all the basic rays. 



op»«»p««o«op»««n 

• •a • • • a ci • 3 «f mcf mtf • 

• 09 • D aro 09 » • pm • d • 

a • • • CK>«ci • kef mm mcf o 
• /» D >• or) • • • p^po • D p 
cf o • ci«r kef •AO. • Q.* • 

(a) 

(b) 



FIG. 13: Examples of ground-state structures at the bound- 
ary between phases 4 and 9: (a) m/ms — 1/5, (b) m/nis = 
1/4. Unit cells are indicated. 



t<mc/i>c/n^i(o^t^t/9^i(mc/iidm 



FIG. 14: An example of a ground-state structure at the 
boundary between phases 4 and 6. Lines separate domains 
of phases 4 and 6. 



VI. GROUND-STATE PHASE DIAGRAMS AND 
FRACTIONAL MAGNETIZATION PLATEAUS 

To make our results more usable, we present all the 
types of ground-state phase diagrams in the plane {h, J2). 
There are three types of diagrams if Ji < (Fig. 18) 
and four types if Ji > (Fig. 19). In these diagrams. 



the lines between neighboring regions correspond to the 
3-faces and the points where three or more regions con- 
verge correspond to the 2-faces. If there is a first-order 
phase transition between the neighboring phases, then 
the line separating the corresponding regions is solid red; 
if the transition is continuous, then the line is solid black. 
The dash-dotted green line between regions 4 and 9 corre- 
sponds to a jump along with a continuous transition. Re- 
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TABLE IV; Basic rays and ground-state configurations for the 2-faces of the full-dimensional regions for the Ising model on 
the extended SS lattice {h ^ 0). 





Ground-state configurations 




Full-dimensional 


Coordinates 




Conditions for existence 


2- face 


for the 2- face 


Disorder 


structures 


in the plane {h, J2) 


in the plane (ft, J2) 


{ri,r2} 


op om 90 mm , mo om n oo o* 

cfo tfo cf* 0^ \ domm II oo wo 





1, 1, 3 


ft = 0, Jl = 




J2 < 0, J3 < 


{ri,r4} 


op mp mm i op n oo om mm 
cfo *o *• 1 *o cf» II oo #0 •• 





1, 1, 9 


(0,-2Ji) 




Jl < 0, J3 < 


{ri.r?} 


om mp mm om mm n oo om mm 

«0 Cf* 0^ II oo AO 


1 


1, 3, 10 


(4Ji,0) 




Jl > 0, J3 < 


1^1,^8} 


mp mm om mm , op om ,, oo o* 

0« 1 «0 VP II OO AO 


1 


1, 9, 10 


ft = J2, Jl = 




J2 > 0, J3 < 


{ri,rio} 


mp am mp mm am n oo o* 
•fo •o cJ» 0» II oo wo •• 


1 


3, 9, 10 


(2Ji,2Ji) 




J3 < — Jl < 


{r2,r3} 


op mp om mp om mm oo mo mm 
cfe do *o d» *• *• II oo #0 •• 


1 


1, 1, 8 


ft = 0, J3 = - Jl 




Jl < 0, J2 < 


{r2,r6} 


op mp om mp om mm mo mm mm 

cfo 0o «o o« «r* II 00 o* 


2 


1, 7, 8 


ft = 2J3, Jl = 




J2 < 0, J3 > 


{r2,r7} 


o* mp om mm 
•roai«««r* 


2 


1, 3, 7 


ft = 4Ji, J3 = 




Jl > 0, J2 < 


{r2,rg} 


op mp om mp om mm n mo om mm 

cfo do *o d» *• *• II wo #0 o» 


2 


3, 7, 8 


ft = 2Ji, J3 = Jl 




Jl > 0, J2 < 


{r3,r4} 


op op mp mm mm n oo mo mm 

do tfo «o d* II oo 00 


2 


1, 1, 5 


(0,-2Ji -2J3) 




< J3 < -Jl 


{r3,r6} 


•o 
•o 


2 


1, 5, 8 


(2Ji -|-2J3,0) 




< -Jl < J3 


{r4,r8} 


op mp mm mm 
«ro«od««« 


2 


1, 5, 9 


h = J2 + 2J1, Jz = 





-J2 < 2Ji < 0, J2 > 


{r5,rs} 


op mp om 
*o •o *o 


2 


4, 5, 9 


ft = J2 -2Ji, J3 = 





< 2Ji < J2, J2 > 


f T 


op mp om mp mm n mo om mm 

iTo Co d* di II AO MO o« 


2 


3, 4, 5 


(2J3,2Ji -2J3) 




< J3 < Jl 


{rs.rio} 


mp om mp mm^ n oo o* mm^, a 
«fo •o d» d» ||oo»o»»o» 


2 


3, 4, 9 


(-2J3,2Ji) 




-Jl < J3 < 


{r6,r7} 


om mp mm om mm n mo mm mm 

i(o di d* i(* «• II vo o* 


2 


1, 6, 7 


{at I T A\ 

(4Ji -1- 2J3,0) 




Jl > 0, J3 > 




op mp om mm om mm n mo mm mm 

iTo «rp «o d* «r« II vo o« 


2 


1, 5, 6 


ft = J2 + 2J3, Jl = 





J2 > 0, J3 > 


{re, rg} 


•o 

#o o# 


2 


5, 6, 7, 8 


(2J3,0) 




< Jl < J3 


{r7,r8} 




2 


1, 4, 6, 10 


ft = J2 + 4 Jl , J3 = 





Jl > 0, J2 > 


{r7,r9} 


om mp mm om mm n om mm i mo 

di d* II AO o* \ mo 


2 


3, 4, 6, 7 


(4Ji ^2J3,0) 




< J3 < Jl 


{r7,rio} 


om mp mm^ o** n om mm i oo 
•p di d* ||«o**o« loo 


1 


3, 4, 10 


(4Ji +2J3,-2J3) 




-Jl < J3 < 




op mp om mm om mm n mo om mm 

ifo •o *o d» *• II wo #0 o» 


2 


4, 5, 6 


ft = J2 4- 2 Jl , Jl = 


J3 


Jl > 0, J2 > 


{r8,rio} 


mp gm mm^ gm^f ,, oo om mm mm^^ 
iTo Co d* mm ||oo«o««o* 


2 


4, 9, 10 


ft = J2, Jl ~ — J3 




Jl > 0, J2 > 


{r9,r^} 


•o o« 
•o mo 


2 


3, 5, 8 


Jl = J3, J2 = 




Jl > 



"Configurations marked by asterisk enter the structures in blocks shown in Fig. 4. 
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FIG. 15: An example of a ground-state structure at the 
boundary between phases 1, 1, and 5 (2-face {r3,r4}), which 
shows that the disorder at this boundary is two-dimensional. 






FIG. 16: Example of a ground-state structure at the boundary 
between phases 3, 4, and 9 (2-face {r5,rio}). The rectangles 
formed by the ferromagnetic chains together with "open" and 
"solid" squares are shown. 
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FIG. 17: Examples of ground-state structures at the bound- 
ary between phases 4, 9, and 10 (2-face {rg, rio}), which show 
that the disorder at this boundary is two-dimensional. In each 
pair of gray circles, one circle should be open and another one 
— solid. Both structures can mix with each other. 



gions 6, 7, and 10, which give rise to an 1/2-plateau, are 
colored and region 4, which gives rise to a 1/3-plateau, 
is shaded. 

As one can see from Fig. 19 and Table II, the width of 
the 1/2-plateau generated by phase 10 is equal to 2 J2 and 
4Ji for sequences of phases 3-10-1 and 9-10-1; for 
all the rest of sequences that give rise to this plateau its 
width is equal to |4J3|. The widths of zero plateaus and 
fractional ones for various sequences of phases as well as 
conditions for their existence are presented in Table V. 

Which sequence of phases corresponds to the magne- 
tization process in ErB4? A single 1/2-plateau has been 
observed in this compound, therefore, only five initial 
sequences in Table V are possible. However, since in- 
teractions Ji and J2 are antiferromagnetic and approx- 
imately equal and interaction J3 is ferromagnetic and 
relatively large, ^ only the sequence 3-10-1 remains. 
This contradicts the statement of Ref. about the mag- 
netic structure of ErB4 for zero field. According to this 
reference, there should be structure AFl rather than the 
Neel phase (AF3). The structure AFl, in the ground- 
state phase diagram (Ji = J2, h = 0) presented in this 
reference, should correspond to our structure 5, however, 
some other structure is depicted there; maybe this is a 
simple inadvertence. 

On the basis of the ground-state structures that we 
have found previously the authors of Ref. [l^ have con- 
structed a ground-state phase diagram for Ji = J2 > 



and for arbitrary J3 and then obtained numerically a 
spin supersolid phase in the quantum model with strong 
Ising anisotropy. They argue that this ground state ex- 
ists in ErB4 which, without magnetic field, has magnetic 
structure 5. However, structure 5 is possible only under 
the condition J3 > Ji > which looks unrealistic. It 
seems that experimental results in favor of this structure 
in ErB4 without magnetic field^'^^ and the experimental 
results yielding that interaction J3 is ferromagnetic and 
largeS. are contradictory. 

The structure with m/nis = 1/2, presented in Ref. 
for TniB4 is erroneous as well. It is a mixture of struc- 
tures 1 and 5. In the model under consideration, this 
structure exists at the boundary between phases 1 and 5 
where Ji < 0. The authors of Refs. Handle state that, in 
TmB4, interaction Ji « J2 is antiferromagnetic, there- 
fore, it is unlikely that this structure could occur in this 
compound. 

It is also unlikely that, in TmB4, the structure (3,7)4 
with m/nis = 1/9 presented in Ref. iJi (Fig. 20) could 
exist. In comparison with the structure (3,4)4 that we 
proposed in Ref. 0, all the antiferromagnetic chains in 
this structure are shifted. In the model under considera- 
tion, the structure (3,7)4 exists at the boundary between 
phases 3 and 7, that is, for positive values of Ji and J3 
but negative values of J2. 

In Ref. d the following experimental values of interac- 
tion parameters for TmB4 are presented: Ji = J2 = 
0.85 K and J3 = 0.3 K. These correspond to the 
sequence of phases 3-4-6-1, however, then the 
1/3-plateau is approximately twice as wide as the 1/2- 
plateau. [If J3 were ferromagnetic (J3 = —0.3 K), then 
the 1 /2-plateau would be twice as wide] . The 1 /3-plateau 
has not been observed in TmB4. Maybe the longer-range 
interactions, giving rise to the fractional plateaus 1/6, 
1/7..., remove at the same time the 1/3-plateau. 

In Ref. [13, to explain the appearance of fractional mag- 
netization plateaus in TmB4, a spin-electron model has 
been considered. For the 1/2-plateau the structure 6b 
(or symmetric one) has been obtained. The SS diagonals 
are not shown in the figures of Ref. [13, therefore we do 
not know whether the structures obtained there are of 
the same type as in Ref. 4 or as in Ref. 9 (Fig. 20). It 
is easy to show that the interaction J4 (the next-nearest 
neighbors along edges) lifts the degeneracy in phase 6 
and stabilizes one of the structures 6a or 66. The fer- 
romagnetic (antiferromagnetic) interaction J4 stabilizes 
the structure 6oi2iii (66). 

It is interesting that in the case where all the three 
interactions are antiferromagnetic (Ji > 0, J2 > 0, and 
J3 > 0) the full-dimensional phases for h ^ (1, 3, 4, 5, 
and 6) can be considered as one-dimensional ordering of 
ferro- and antiferromagnetic chains. A question arises: In 
order to study the effect of longer-range interactions on 
ground states, is it possible to consider a one-dimensional 
model with effective interactions between chains instead 
of the two-dimensional one?— 
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FIG. 18; Ground-state phase diagrams of the Ising model on the extended SS lattice for Ji < 0. Black and red lines correspond 
to continuous phase transitions and jumps, respectively. (See also Table and Figs. 5 and 6). 




FIG. 19: Ground-state phase diagrams of the Ising model on the extended SS lattice for Ji > 0. Black and red lines correspond 
to continuous phase transitions and jumps, respectively; a green dash-dotted line denotes a jump together with a continuous 
phase transition. The regions which give rise to an 1/2-plateau are colored. The region of 1/3-plateau phase is shaded. (See 
also Table and Figs. 5 and 6). 



VII. SUMMARY 



To conclude, we have determined a complete solution 
of the ground-state problem for an Ising model on the ex- 
tended SS lattice with an additional interaction J3 along 
the diagonals of "empty" squares. We have used the 



method of basic vectors and basic sets of cluster config- 
urations that was proposed in our previous works. Here, 
however, we generalize the method and consider config- 
urations of two clusters at once. We have constructed 
the ground-state phase diagrams and studied the ground- 
state structures in the full-dimensional regions as well as 
at their three- and two-dimensional boundaries. This 
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TABLE V: The widths of the fractional plateaus for various ways of transitions from the zero-field phase to the ferromagnetic 
one. 



Sequence 


Width 


Width 


Width 




of phases 


of zero plateau 


of 1/3- plateau 


of 1/2-plateau 


Conditions for existence 


3-10-1 


4Ji - Ja 





2J2 


Ji > 0, J3 < 0, < J2 < min[2Ji, -2J3] 


3-7-1 


4Ji - 2J3 





4J3 


< J3 < Ji 


5-6-1 


J2 + 2 J3 





4J3 


< Ji < J3 


8-7-1 


2J3 





4J3 


< Ji < J3 


9-10-1 


J2 





4Ji 


J-i < -Ji < 0, J2 > 2Ji 


3-4-10-1 


4Ji - 2 J2 - 2J3 


3J2 + 6J3 


-4J3 


-Ji < J3 < 0, -2J3 < J2 < 2Ji 


9-4-10-1 


-2 Ji + J2 - 2J3 


6(Ji + J3) 


-4J3 


-Ji < J3 < 0, J2 > 2Ji 


3-4-6-1 


4Ji - 2 J2 - 2J3 


3J2 


4J3 


< J3 < Ji, < J2 < 2(Ji - J3) 


5-4-6-1 


-2 Ji + J2 + 4J3 


6( Ji - J3) 


4J3 


< J3 < Ji, J2 > 2(Ji - J3) 




FIG. 20: Structure (3,4)4 for the 1/9-plateau in TmB4, pro- 
posed in Ref.Q, and structure (3,7)4, proposed for this plateau 
in Ref. Qj. These two structures differ by the positions of the 
antiferromagnetic chains. Unit cells are indicated. 



made it possible to establish that the additional inter- 
action J3 gives rise to an 1/2-plateau. This plateau can 
correspond to three different phases (depending on the 
interaction parameters), two of which are partially disor- 
dered. In addition to the 1/2-plateau, another fractional 
plateau is possible — with the magnetization 1/3. As we 
have shown earlier," this is a single fractional plateau 
in the Ising model on the conventional SS lattice. For 
some relations between the interaction parameters, the 
1/2-plateau can be a single fractional plateau in the Ising 
model on the extended SS lattice as well. Hence, it might 
be reasonable to believe that we have explained the ori- 
gin of a single fractional plateau — with m/m^ = 1/2 — in 
ErB4. As to TmB4, where no 1/3-plateau but an 1/2- 
plateau and a sequence of other fractional plateaus were 
observed, the theoretical explanation of the magnetiza- 
tion curve in this compound requires to study the effect of 
longer-range interactions. However, the knowledge of the 
ground-states at the boundaries of the full-dimensional 



regions of the model under consideration makes it pos- 
sible to draw some conclusions about the effect of such 
interactions. 

The advantages of the analysis of the ground states 
at the boundaries of full-dimensional regions become ob- 
vious when the results of this paper are compared to 
the results of the previous one.^ Having constructed the 
ground-state phase diagram for the Ising model on the 
conventional SS lattice and having analyzed which of 
the ground states at boundaries of full-dimensional re- 
gions become full-dimensional if a small additional inter- 
action J3 is switched on, we have obtained all the full- 
dimensional ground-state structures of the Ising model 
on the extended SS lattice except for structure 8 (since 
it occurs under the condition J3 > | Ji | ) . 

The results obtained here and in Ref. [3 can be useful for 
the numerical studies of the origin of the fractional mag- 
netization plateaus in SS magnets (with Ising anisotropy) 
for nonzero temperaturesji^ 

In the present paper, a development of analytical meth- 
ods for the determination of ground states of lattice-gas 
models and equivalent spin models is presented. These 
methods can be used to investigate the structure of sub- 
stitutional alloys, that is a very important problem (a 
recent initiative of US President Barack Obama-- gives 
evidence of this). In our future paper we will develop the 
method and show how to study the effect of longer-range 
interactions. 
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